Search for Superfluidity in Solid Hydrogen 
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A torsional oscillator study of solid para-hydrogen has been carried out down to 20 mK in search 
for evidence of superfluidity. This work was inspired by the observation of the supersolid phase in 
solid 4 He. We found evidence of a possible phase transition, marked by an abrupt increase in the 
resonant period of oscillation and onset of extremely long relaxation times as the temperature was 
raised above 60 mK. The change in the period for para-hydrogen, in contrast to solid 4 He, is not 
related to superfluidity as it fails a crucial test showing that it is not a consequence of irrotational 
flow. The long relaxation times observed suggest the effect is related to the motion of residual 
ortho-hydrogen molecules in the solid. 

PACS numbers: 66.35.+a, 66.30.-h, 67.80.-s, 67.80. Mg 
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A supersolid is a solid that possesses an element of su- 
perfluidity. The existence of supersolid 4 He was first sug- 
gested 67 years ago 0, while possible underlying mecha- 
nisms were discussed 30 years later Q . Experimental ev- 
idence of supersolid 4 He was reported by Kim and Chan 
(KC) only two years ago in a series of torsional oscillator 
experiments on bulk solid helium 0-0 and solid helium 
confined in porous media 0, . KC demonstrated that 
the torsional oscillator is an ideal technique for the detec- 
tion of solid 4 He acquiring non-classical rotatonal inertia 
(NCRI) 0. NCRI is indicated by a drop in the resonant 
period of the oscillator below the transition temperature. 
The decrement in the period is proportional to the su- 
persolid fraction in the limit of low oscillation speed Q . 
These observations have since been replicated in three 
other laboratories ||. This prompted us to address the 
question of whether or not the supersolid phase is present 
in any system other than 4 He. 

Helium is the most quantum mechanical solid of all 
the elements. This is due to the combination of large 
zero-point energy of individual atoms and weak attrac- 
tive interactions among them. A measure of the quantum 
nature of a substance in the condensed phases, i.e. liq- 
uid or solid, is given by the de Boer parameter, A = 
h/ay/me, where a and e are fitting parameters for the 
inter-particle potential energy Only helium and hy- 
drogen isotopes have A greater than unity (A^He =3.0, 
A iHe = 2.6, A H2 = 1.7, A HD = 1.4, A D2 = 1.2). There- 
fore H2, after 4 He, is the most likely system to exhibit 
superfluidity. Indeed the zero-point motion of solid H2 
at 5.4 K under zero pressure is 18% of the nearest neigh- 
bor distance while in low density solid 4 He it is 26% 
11]. Both of these values are significantly higher than 
the Lindcmann (melting) criterion of approximately 10% 
for classical solids. The fact that the supersolid phase 
in 4 He persists up to at least 135 bar [4|, where it is 
expected to be less quantum mechanical than at lower 
pressures, was further incentive to investigate solid H2. 

We have carried out a series of torsional oscillator mea- 
surements on solid hydrogen. Our experimental configu- 



ration is similar to that used by KC J3j- A sample cell is 
suspended from the cooling stage of a 3 He- 4 He dilution 
refrigerator by an annealed beryllium-copper torsion rod. 
The rod provides a restoring force to keep the cell in tor- 
sional motion and doubles as a filling line to the cell. The 
H2 space in the cell is comprised of an open annulus and 
a central filling line, which are connected by three chan- 
nels that span the diameter of the cell. The annulus has 
an average radius, width, and height of 6.5 mm, 2.3 mm, 
and 4.8 mm, respectively. The radial channels and filling 
line have a diameter of 1 mm, a combined volume that 
constitutes 7% of the total cell volume, and contribute 
about 2% of the rotational inertia of all the H2 in the 
cell. Electrodes on the exterior of the torsion cell are 
used to capacitively drive and detect oscillatory motion. 
The resonant period (to) of the oscillator, is determined 
by the moment of inertia (I) of the torsion cell and the 
torsion spring constant (ft) of the rod, such that tq — 
2-k^/TJk. The mechanical quality factor (Q) of the os- 
cillator at low temperature is 1.1 x 10 6 , which allows for 
the detection of changes in tq of one part in 10 7 . 

Solid H2 samples were grown at saturated vapor pres- 
sure from high purity H2 gas containing very low levels of 
isotopic impurities: less than 10 ppm of both hydrogen 
deuteride (HD) and deuterium (D2). However, for pure 
H2 crystals there exist impurities of a second kind. The 
two lowest rotational energy states of each molecule are 
the para- (P-H2) and ortho- (0-H2) states, having angular 
momentum J = and J = 1, respectively. The equilib- 
rium concentrations of P-H2 (1 - x) and 0-H2 (x) at 300 K 
are respectively 0.25 and 0.75. At temperatures (T) be- 
low 4 K the equilibrium value of x is approximately zero 
and the solid is essentially pure P-H2. However, the con- 
version rate from 0-H2 to P-H2 in the solid is very slow 
for T > 500 mK (eg. about three weeks of conversion 
are required for x to drop from 0.01 to 0.009). This rate 
can be increased substantially by keeping H2 in the liquid 
phase (with T > 13.8 K, the triple point temperature) in 
the presence of a high surface area, magnetic material. In 
this way we attain x < 0.01. This process is irreversible 
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so that the average x of each sample can only decrease 
with time. Concentrations are determined by keeping 
the mixing chamber at 20 mK and measuring the tem- 
perature difference between the torsion cell and mixing 
chamber. Based on the temperature difference and the 
thermal conductance of the Be-Cu torsion rod, we can 
infer the amount of o-p conversion heat being released 
and hence x within the H2 sample. The smallest mea- 
surable temperature gradient corresponds to x = 0.01. 
All H2 samples discussed in this work have x < 0.01 (i.e. 
no observed gradient). To limit the uncertainty in the 
absolute value of x, we measured the thermal conductiv- 
ity of our Be-Cu torsion rod and found good agreement 
with extrapolated values from high temperature ( T > 1 
K) data previously reported 0] . 

In Fig. \Ql we have plotted the resonant period as a 
function of temperature for several samples. The empty 
cell period at 1 K (tik) is 709,700 ns, which decreases 
by only 3 ns, or five parts in 10 6 , as the cell is cooled to 
20 mK. This change is determined by the T-dependence, 
if any, of / and k. Based on earlier work the small 
change in tq depicted in Fig. is due to small changes 
in / and/or n rather than large, nearly compensating 
changes of each. For the H2 sample in Fig. ^ the pe- 
riod increases from the empty cell value by 2682 ns upon 
solidification, indicating that the cell is 88% full. The 
mass loading due to the heavier and denser isotope, HD, 
is 4675 ns, corresponding to 91% filling. Each data set 
has been shifted vertically for easy comparison. The cor- 
responding values of t\k for each are noted in the figure. 
There is a slight difference in the slope of each curve 
shown. This modification of the "background" when the 
sample is introduced into the cell is always observed in 
torsional oscillator studies of liquid helium films , as 
well as in the solid helium work 0, HJ • 

The most striking feature in the data of Fig. Q]i is the 
clear difference between H2 and HD. The abrupt drop in 
To below 180 mK for the cell containing H2 is absent for 
HD. This discrepancy is similar to that found between 
4 He and 3 He, and indicates that the drop in the period 
could be a signature of NCRI, as it is seen only in bosonic 
solids. The net change in To is obtained by subtracting 
the measured period from values extrapolated from high 
temperature, as indicated in the figure. While the pe- 
riod drop is more than an order of magnitude smaller 
than that found in solid 4 He in a similar geometry U its 
temperature dependence is similar [Fig.^>]. The change 
in To first becomes measurable below 180 mK, and in- 
creases first gradually with decreasing T and then much 
more rapidly prior to saturation near 60 mK. This simi- 
larity in the T-dependence also suggests that the period 
drop in H 2 could be a signature of NCRI. 

However, there are four puzzling features in the H2 
data. First, for solid 4 He the drop in the period is ac- 
companied by a minimum in the amplitude of oscillation 
at the temperature where tq changes most rapidly 



There is no such minimum in oscillation amplitude near 
the transition region in the H2 experiment. Second, for 
4 He the measured NCRIF decreases once the maximum 
oscillation speed of the sample exceeds 10 /zm/s 0- This 
speed has been interpreted as the critical velocity of su- 
perflow. There is no evidence of a critical velocity in solid 
H 2 . We found the period drop to be independent of os- 
cillation speed up to 500 /im/s, the maximum within the 
linear response range of the oscillator. 

Third, the apparent supersolid onset temperature (180 
mK) and saturation temperature (60 mK) are remark- 
ably close to the 230 mK and 50 mK values found for 
solid 4 He The similarity in these characteristic tem- 
peratures is inconsistent with the substantial difference 
in K^He and A#2- For solid 4 He in Vycor glass, very low 
3 He impurity levels can significantly enhance the onset 
temperature |5j. While the Tq for a sample containing 
0.2 ppm of 3 He is about 175 mK, a sample with 10 ppm 
has an onset temperature around 370 mK. This effect has 
also been observed in bulk 4 He , thus a small amount 
of HD impurities may alter Tc for solid H2. However, 
we do not know the exact HD impurity level in our solid 
H2 samples. The reason is that HD is heavier than H2 
and hence preferentially adsorbs onto surfaces, which are 
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"Open Annulus Cell": 

—X— Empty, t m = 709,700 ns 

— HD, t ik = 714,375 ns 

— ♦— H2[.r< 1%], r IK = 712,382ns 

"Blocked Annulus Cell": 

— A— H2 [ v < 1%], r IK = 783,629 ns 












































"Open Annulus Cell": 

— ♦— H2 [88% full,.r< 1%] 

"Blocked Annulus Cell": 

H2 [80% full, x < 1%] 
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FIG. 1: a.) T-dependence of the resonant period for both 
the "open" and "blocked" cells. All of the data shown were 
obtained at maximum rim velocities less than 40 fim/s. For 
the open annulus, data is plotted for the empty cell, cell con- 
taining HD, and cell containing H2. H2 data taken using the 
blocked cell exhibits the same qualitative features, including 
long relaxation times, b.) T-dependence of the absolute pe- 
riod change from the high- T background, normalized by the 
total H2 mass loading ("NCRIF," the apparent NCRI frac- 
tion). The T-dependence is similar to that observed in solid 
4 He. However, the same effect is seen in the blocked cell. 
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abundant in the H2 gas handling system leading into the 
torsion cell. Specifically, the capillary system contains a 
chamber filled with a high surface area, FeO(OH) com- 
pound. Thus, we expect that the actual HD concentra- 
tion in each sample is much less than the quoted value of 
10 ppm. Since we do not know with certainty the exact 
HD concentration we cannot validate the 180 mK onset 
and 60 mK saturation temperatures of NCRI. 

The final and most unusual feature of the H2 data, 
compared to 4 He, is the exceedingly long relaxation time 
(t) of the period reading [Fig. 0. It is expected that the 
oscillator comes into thermal and mechanical equilibrium 
in the following manner. The e _1 time constant required 
for thermal equilibration increases with decreasing T and 
typically reaches 15 minutes below 50 mK. Since the me- 
chanical Q is on the order of 10 6 , the necessary time 
to allow for mechanical equilibrium is about 10 minutes. 
Therefore it is necessary to wait at least 75 minutes af- 
ter each temperature step before collecting To data. This 
protocol was used in the solid helium experiments to en- 
sure reproducible data upon warming and cooling. How- 
ever, in gathering data on H2 we observe lvalues as long 
as 36 hours just above 60 mK. There is no such evidence 
of relaxation in the HD sample or for the empty cell, for 
which To equilibrates within 75 minutes at all T. We will 
return to the issue of lvalues below. 

We have carried out a crucial control experiment to 
ascertain whether the observed effect in solid H2 is a sig- 
nature of superfluidity. We constructed another sample 
cell, identical to the first except for the addition of a bar- 
rier at one point in the annulus. Its purpose is to make 
the path for the irrotational superflow much more tortu- 
ous, leading to a predictable decrement in NCRI 0,0], 
thus reducing the drop in the period. This control exper- 
iment provides a convincing case for the supersolid phase 
of 4 He. It was found that when the annulus was blocked, 
the drop in To was reduced by a factor of 70 from the 
open annulus value Q. This reduction is in agreement 
with what is expected for irrotational superflow |l6J • The 
result of this control experiment for solid H2 is plotted 
alongside the data of the "open annulus cell" in Fig. QJt. 
The period drop below 180 mK is still present for the 
"blocked annulus cell." In fact, we find the magnitude 
is larger in the latter. One possible reason for this is 
that, while each sample has x < 0.01, the exact ortho- 
concentrations of the two samples may differ. Regardless, 
the result of the control experiment clearly indicates that 
the H2 result is completely different from that found in 
solid 4 He and is not related to superfluidity. Wc note 
that, as in the open cell, long relaxation times are ob- 
served just above the saturation temperature of 60 mK. 

Long equilibration times are typical in studies of 0-H2 
diffusion in solid P-H2 |l7^ . In the dilute limit where x 
is less than several percent, isolated o-H 2 molecules (sin- 
gles) cluster together below 1 K even though thermally 
activated diffusion is severely limited. The 0-H2 impuri- 



ties effectively propagate throughout the lattice via the 
transferral of J from molecule to molecule [l8j |. rather 
than by particle exchange. The "equilibrium" state of 
the solid thus involves a temperature dependent distri- 
bution of singles, pairs (two neighboring 0-H2 molecules) 
and larger clusters of 0-H2. At several kelvin, all of the 
0-H2 molecules are randomly distributed as singles. As 
the temperature is lowered the concentration of isolated 
singles decreases in favor of the formation of pairs and 
larger clusters. Since 0-H2 and P-H2 have different den- 
sities, in the limit of zero temperature the system prefers 
to be in a macroscopically phase-separated state. 

Clustering was first observed by monitoring the time 
evolution of the nuclear magnetic resonance (nmr) ab- 
sorption spectra for 0-H2 singles and pairs after rapid 
warming or cooling of solid H2 crystals |lj| . Representa- 
tive growth (above 1 K) and decay (below 1 K) rates of 
the singles spectrum [2(j are presented in Fig.|2] The sin- 
gles decay rate exhibits a non-monotonic T-dependence, 
decreasing dramatically between 300 mK and 100 mK. 
Well below 100 mK the decay of singles is on a time 
scale of minutes rather than hours, rendering the mea- 
surement impossible. The pair spectrum also shows an 
enhancement in its decay rate at low T but time con- 
stants are at least one order of magnitude longer, eg. on 
the order of 10 hours at 25 mK [2lJ. Also included in 
Fig. are relaxation times obtained from simultaneous 
thermal conductivity measurements |2(i| . The time evo- 
lution of the thermal conductivity depends strongly on 
T. For T > 300 mK the time dependence is described 
by a simple exponential. However, contrary to the sin- 
gles decay rate, a sum of two exponentials best fits the 
conductivity data below 300 mK. This difference may be 
related to the presence of 0-H2 clusters at lower temper- 
ature. Since it is difficult to resolve the spectra of triples 
or larger clusters j^, nmr studies cannot quantify the 
degree of mobility of larger 0-H2 clusters that form at low 
T. Their presence, however, has been detected for solid 
H2 crystals held at 25 mK for only 24 hours |23|. 

Our torsional oscillator data are obtained using the fol- 
lowing protocol. Upon cooling H2 samples to 20 mK, one 
to two weeks are allowed for equilibration, during which 
time To drops smoothly until finally stabilizing to within 
0.05 ns, our limiting resolution. We note that the rate of 
irreversible o-p conversion is enhanced for T < 200 mK 
since local a;-values are higher in clustered regions |24j . 
However, upon stabilizing at 20 mK, we observe no long 
term drift in the To-reading. After complete equilibration 
the temperature is raised in successive steps, for each of 
which t is measured as a function of time [Fig. |21 inset] . 
We find for all samples that t equilibrates quickly for 
any T < 60 mK. However, further increase in T leads 
to extremely slow relaxation. Near and above 180 mK 
we find that t shortens, returning to a value of about 
15 minutes, the usual e _1 equilibration time. When the 
temperature sweep is complete the system is returned 
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nmr [0.4 - 3% o-H2] 
th cond [0.4 - 3% 0-H2J 
tors osc [0.5% o-H2] 
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FIG. 2: T-dependence of measured relaxation times, includ- 
ing data from the literature |20|. The sharp increase in t 
coincides with the observed period shift. At T ~ 180 mK, a 
minimum in t is observed for both the torsional oscillator and 
thermal conductivity measurements. Inset: Time dependence 
of ro at several different temperatures. The curves have been 
shifted vertically for clarity. The longest relaxation times are 
observed at T ~ 80 mK. 



to 20 mK and allowed to re-equilibrate. Our measured 
i-values are presented in Fig. [21 

Associating the motion of 0-H2 impurities with the ob- 
served time dependence of ro in the present work explains 
the absence of this unusually long relaxation in HD, 4 He, 
and He. In our experiment we allow clustering to en- 
sue for up to two weeks at 20 mK. Based on measured 
decay times in the literature, no isolated 0-H2 singles or 
pairs remain in our samples at the time we begin each 
temperature sweep. Thus, we are able to follow the un- 
clustering process upon gradual warming of the torsion 
cell. 

The moment of inertia and hence To are extremely sen- 
sitive to the radial density profile of the sample. Changes 
in the period require the rearrangement of 0-H2 within 
the torsion cell in such a way as to reduce or enhance 
the density at large radii. We think the most likely ex- 
planation for the observed drop in ro below 180 mK is 
the formation of large clusters of nearly pure 0-H2 in the 
center of the cell (fill line and radial channels), increasing 
the purity of P-H2 in the annulus. Since the density of 
the 0-H2 is 1.7% higher than that of p-Eb, the clustering 
of 0-H2 toward the center has the net effect of reducing 
7. An extreme case is where all of the 0-H2 molecules 
cluster to the center of the torsion cell (the filling line). 
In this situation, the reduction of To is 1.7 x 10~ 4 of the 
total H2 mass loading. This is on the same order of our 
observations [Fig. H], i.e. 5.5 x 10 -4 and 9.0 x 10~ 4 in 
the open and blocked cells, respectively. It is known that 
solid H 2 does not wet many metal substrates |25j . If the 



contact angle of the sample to the cell walls is different 
for the mixed and phase-separated phases (eg. P-H2), I 
will show a change upon phase separation. This is one 
possible explanation for larger than expected change in 
To that we observe. 

In conclusion, a resonant period shift in our torsional 
oscillator containing solid H2 is observed. While the phe- 
nomenon shares some features with the supersolid tran- 
sition in 4 He, there are several dissimilarities. Most no- 
table is the presence of the period drop in both the open 
and blocked cells, proving that the effect is not related to 
superfluidity. We believe the abrupt rise in ro accompa- 
nied by a sudden increase in t near 60 mK is consistent 
with a transition from a phase-separated configuration to 
one that is un-clustered at high temperature. 

We thank E. Kim and W. N. Hardy for their advice. 
Special thanks to H. Meyer for many enlightening discus- 
sions. Support comes from NSF Grant DMR 0207071. 
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